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ABSTRACT 


We report on a multi-wavelength observational campaign which followed the 
Earth’s transit on the Sun as seen from Jupiter on 5 Jan the 2014. Simultaneous 
observations of Jupiter’s moons Europa and Ganymede obtained with HARPS from 
La Silla, Chile, and HARPS-N from La Palma, Canary Islands, were performed to 
measure the Rossiter-McLaughlin effect due to the Earth’ s passage using the same 
technique successfully adopted for the 2012 Venus Transit (Molaro et al. 2013). The 
expected modulation in radial velocities was of « 20 cm but an anomalous drift 
as large as « 38 m i.e. more than two orders of magnitude higher and oppo¬ 
site in sign, was detected instead. The consistent behaviour of the two spectrographs 
rules out instrumental origin of the radial velocity drift and BiSON observations rule 
out the possible dependence on the Sun’s magnetic activity. We suggest that this 
anomaly is produced by the Opposition Surge on Europa’s icy surface, which ampli¬ 
fies the intensity of the solar radiation from a portion of the solar surface centered 
around the crossing Earth which can then be observed as a a sort of inverse Rossiter- 
McLaughling effect, in fact, a simplified model of this effect can explain in detail most 
features of the observed radial velocity anomalies, namely the extensions before and 
after the transit, the small differences between the two observatories and the presence 
of a secondary peak closer to Earth passage. This phenomenon, observed here for the 
first time, should be observed every time similar Earth alignments occur with rocky 
bodies without atmospheres. We predict it should be observed again during the next 
conjunction of Earth and Jupiter in 2026. 
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1 INTRODUCTION 

Transits of Venus and Mercury in front of the Sun are major 
historical events but also other transits can be seen in the 
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solar system from other planets each time the heliocentric 
conjunctions take place near one of the nodes of their orbits, 
with the exception of the innermost Mercury. In particular, 
the Earth can be seen transiting in front of the Sun from 
other planets. These are rare events which were predicted in 
detail by |Meeu's| ( |1989[ ). For instance, the Earth will be seen 
transiting the Sun from Mars only in 2084. As seen from 
Jupiter, a transit took place on 5 January 2014. Next passage 
will be grazing and will occur in 2026. During these transits, 
the integrated solar light can be recorded as it is reflected 
by the planets from which the Earth is seen transiting in 
front of the Sun, offering a surrogate direct watch as we 









2 P. Molaro, M. Barbieri, L. Monaco, S. Zaggia and C. Lovis 


showed with the observation of the Venus Transit of 6 June 
2012 when we followed the transit as if it were seen from the 
Moon ( Molaro et al.|20T3 l. 

We planned an observational campaign to observe the 
Earth’s passage in front of the Sun that took place in Jan¬ 
uary 2014. One of the motivations for this observational 
campaign was the detection of the Rossiter-McLaughlin 
(RM) effecton which we report in this work. The RM is a 
radial velocity drift caused by the distortion of the stellar 
line profiles due to the occultation of the rotating stellar disk 
by an intervening body. The effect was first predicted by 


Holt (1893), and discovered by Schlesinger (19111, and later 
confirmed by Rossiter (19241 and McLaughlin ( 1924[) in the 
eclipsing binaries /3 Lyrae and Algol, respectively. |Schnei-| 


der (20001 suggested that the transit of a planet could also 


be detected in the line profile of high signal-to-noise stel¬ 
lar spectra of rotating stars, and a Jupiter-like planet was 


first observed in HD 209458 by Queloz et al. 120001 with 
an amplitude of ± 30 m s“^. The detection of the RM ef¬ 
fect provides information on the planet radius, the angle A 
between the sky projections of the orbital axis and the stel¬ 
lar rotational axis. Since then about 90 other Jupiters have 


been observed, often with very tilted orbits (Fabrycky & 
Winn|[^09[ |Triaud et al.|20l0| [Brown et al.||2012| [Albrecht 


2012). The smallest RM effect detected is due to the Venus 


Transit in front of the Sun of 6 June 2012 by [Molaro et al.| 
(2013) who used the integrated sunlight as reflected by the 


Moon at night-time to record about half transit by means 
of the high precision HARPS spectrograph at the 3.6m La 
Silla ESO telescope. The observations performed in corre¬ 
spondence of the passage of Venus in front of the receding 
solar hemisphere showed that the planet eclipse of the solar 
disk was able to produce a modulation in the radial velocity 
with an amplitude of « -1 m The radial velocity change 
is comparable to the solar jitter and is more than one order 
of magnitude smaller than that of extra-solar hot Jupiters. 

The amplitude of the radial velocity anomaly stemming 
from the transit is strongly dependent on the projected ra¬ 
dius of the eclipsing body and on the component of the star’s 


rotational velocity along the line of sight ( 

Ohta, Taruya & 

[Suto[[2005[ [Gimenez|[2006[ [Gaudi & Winn| 

20071). A transit 


across a star with high projected rotational velocity pro¬ 
duces a radial velocity signature larger than across a slow 
rotator. The radial velocity drift AV), is given by: 


AK = sin/. 


( 1 ) 


where is the stellar angular velocity, Sp is the pro¬ 
jected position of the planet on the stellar surface Sp = 
{Xp + Zp)^^^, I, is the inclination between the stellar spin 
and the y-axis and k = Rp/R, isthe ratio between the planet 
and stellar radii ( Ohta, Taruya &: Suto|2005 ). 

During the Earth’s transit of 5 January 2014 the pro¬ 
jected size of the Earth was about 1.3 xlO”"^ of the solar 
disk. Assum ing a solar rotation v elocity of v sin 70 = 1.6 ± 
0.3 km s ^ ( Pavlenko et al.|2012 | the expected RM effect is 
of the order of ± 20 cm s“^. Furthermore, our Moon is also 
transiting the solar surface but with a delay of about four 
hours with respect to the Earth. This type of configuration 
should be quite common in transits of extrasolar planets 
which likely have also their own moons. . The expected RM 
effect due to the Moon’s occultation is of only few cm s“^. 



Figure 1. Composite image of the Sun with Earth and Moon as 
seen from Europa at 19 UT of 5 January 2014. The sizes are in 
scale with the Earth’s size of 4.2 arcsec and the solar disk’s size 
of 369 arcsec. The solar image is from SDO/NASA HMI Inten- 
sitygram at 617.3 nm on 5 Jan 2014 and shows prominent solar 
spots in the approaching solar hemisphere. 
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Figure 2. Front view of the Jovian system on 5 January 2014 
from an observer on the Sun or on the Earth. Prom Jupiter the 
Earth transit started at MJD 56662.70, while its moons arrive 
somewhat in advance to the alignment, of about one hour for 
Ganymede and about 30 minutes for Europa, while lo went behind 
Jupiter during the event. 


2 OBSERVATIONS 
2.1 Timing of the transit 

In Fig. the Earth, the Moon and their trajectories are 
shown as they would appear to an observer on Jupiter (or 
on one of its moons) on 5 January 2014. First contact was 
at solar latitude of -23.8° while the exit was at -35.7°. The 
heliographic latitude of the centre of the disk, the solar Bo 
angle, was of -3.6° and therefore the Sun was showing the 
south pole to Jupiter with an inclination of 6° East of the 
solar axis. From the Jovian system, the black disk of the 
Earth was of 4.2 arcsec while the whole solar disk was of 369 
arcsec. The total duration of the passage was of 9h40m. 

Jupiter itself is not a good sunlight reflector due to its 
high rotational velocity and to the turbulent motions of its 
atmosphere. Its major solid moons are better reflective mir¬ 
rors. The geometrical configuration of the Jovian system is 
illustrated in Fig. [^from an observer on the Sun. The timing 
of the Earth’s transit varies from one moon to another. In 
January the moons were seen approaching Jupiter and there¬ 
fore arrived at the alignment slightly before the planet. The 
Earth transit on the reference frame of the jovian system 
started at MJD 56662.70 from Jupiter, but it was seen by 
about 30 minutes in advance from Europa and about one 
hour from Ganymede. 
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On January 2014 Jupiter could be seen at best from 
the northern hemisphere, but there was not a suitable site 
where Jupiter could have been observed during the entire 
10-hour transit. The moon Europa was the best suitable 
replacement for Jupiter, providing the most extended cov¬ 
erage of the transit for about 6 hours from La Palma and 
offering a limited possibility from La Silla to follow for « 
1 hour the end of the transit. From La Silla it was possi¬ 
ble to observe the beginning of the night when Jupiter was 
rising at 20 degrees over the horizon, but remaining always 
quite low and reaching 35 deg at the end of the Transit. 
The transit could not be observed from Manna Kea either 
and high resolution facilities that could deliver very precise 
radial velocity measurements were not available in other as¬ 
tronomical sites. Thus, La Palma and La Silla were the only 
sites where the phenomenon could be followed with high res¬ 
olution spectrographs suitable to deliver the required radial 
velocity precision. 


2.2 The observations 


The observations comprise a series of spectra taken with 
both HARPS-N and HARPS of the Jupiter’s moons Europa 
and Ganymede covering the range from 380 to 690 nm. At 
the epoch of the observations Europa and Ganymede were 
fully illuminated and had a visual magnitude of 5.35 and 
4.63 mag and apparent diameters of 1.02 and 1.72 arcsec, 
respectively. The integration times of the observations were 
60 or 120 s and delivered a signal-to-noise ratio of « 200 
each at 550 nm with a resolving power of R = A/AA « 
115000. The two spectrographs at La Silla and La Palma 
are twins. Both are in vacuum, thermally isolated, stable 
and equipped with an image scrambler which provides a 
uniform spectrograph-pupil illumination which is essential 
for high precision radial velocity observations. HARPS was 
able to deliver a sequence of observations with a dispersion 
of 0.64 m over a 500-day baseline for the radial velocity 
curve of an extra-solar planetary system composed by three 
Neptune-mass planets ( Lovis et al.||2006 l. 

The observations started as soon as Jupiter’s moons 
became observable from the two sites. We started observ¬ 
ing Ganymede from both telescopes on the night preced¬ 
ing the transit to determine the pre-transit characteristic 
solar radial velocity. At La Palma the observations began 
on 2456661.983 MJD till 56662.265 MJD and La Silla on 
56662.088 MJD till 56662.321 MJD. 


The following night we observed Europa from both 
telescopes to cover the second fraction of the transit as 
much as possible. At La Palma observations started at 
MJD 56662.859 and ended at MJD 56663.265, while at La 
Silla observations were taken in the interval between MJD 
56663.070 and 56663.330. 


In the night following the transit we made observations 
of both Europa and Ganymede to determine the post-transit 
characteristic solar radial velocity only from La Silla. Ob¬ 
servations of Europa were taken from 56664.068 MJD to 
56664.167 MJD, followed by a sequence of observations of 
Ganymede till 56664.327 MJD. 



56662.0 56663 56663.2 

\«D 


Figure 3. Radial velocities measured from the Europa spectra of 
5 of Jairuary with HARPS-N observations. The top panel provides 
the RV corrected for the kinematical motions of the observer and 
of the moons, which are shown in the bottom panel with a black 
line. The red dotted line are the observed radial velocities. Notice 
the 5 m oscillation of the Sun which is responsible for high 
frequency variations of an amplitude of ± 1 m 


3 RADIAL VELOCITIES 


We used HARPS and HARPS-N pipelines to obtain the ra¬ 
dial velocities from the observations. The pipelines return 
a radial velocity value from the cross-correlation of the 
spectrum with a G2 V flux template which is the Fourier 
transform spectrometer (FTS) obtained by Kurucz at the 
McMath-Pierce Solar Telescope at Kitt Peak National Ob¬ 
servatory ( Kurucz et al.|[T984[ ). The FTS solar spectrum is 
calibrated onto the telluric emission lines and is known to 
have an offset in the zero point of the order of 100 m s~^ 


(Kurucz et al. 1984 Molaro & Monai 20121. The pipeline 


returns the radial velocity RVp relative to the solar system 
barycenter by taking the apparent position of the target. 
We thus subtracted the barycentric radial velocity correc¬ 
tion, the BERV, which was recorded in the fits headers to 
compute the proper kinematical corrections. These included 
the motions of the observer relative to Jupiter’s moons at 
the instant when the light received by the observer was re¬ 
flected by the moons, but also the radial velocity compo¬ 
nents of the motion of the moons relative to the Sun at the 


instant the light was emitted by the Sun (Molaro & Gen 
|turi6n|[2011| [Lanza fc Molaro||2015[ ). The sunlight reflected 
by Jupiter’s moons is shifted by the heliocentric radial ve¬ 
locity of the moon with respect to the Sun at the time the 
photons left Jupiter’s moon and were shifted by the compo¬ 
nent of the Earth rotation towards the moon at the time the 
photons reach Earth. The latter is the projection of the as¬ 
teroid motion along the line-of-sight adjusted for aberration, 
and comprises both the radial velocity of the moon and the 
component from the Earth rotation. Thus, the radial veloc¬ 
ity is: 


RV = RVp — BERV — {RVmoon-obs -I- RV^oan-O)- (2) 
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MJD 

Figure 4. Radial velocities measured on 4-6 January 2014. Open black circles are observations of Europa from La Palma while color 
squares are observations of Ganymede (cyan) and Europa (red) from La Silla. A constant offset of 107.5 m as measured far out from 
the transit is taken as the instrumental baseline and is subtracted from the data. The vertical dashed lines mark the expected ingress 
and egress of the Earth’s transit as seen from Europa 


The quantities are computed by using the JPL hori¬ 
zon ephemerides The average rate in the radial velocity 
change of Ganymede and Europa are of about 11 m 
and 12 m per minute, respectively. During the expo¬ 
sure of one or two minutes this velocity change produces 
some spectral smearing. However, we apply the corrections 
to mid-exposure values and since the spectral smearing is 
symmetrical to a good approximation it does not result into 
a net shift in the measured radial velocities. 

Fig-i in the top panel shows the corrected radial ve¬ 
locities for the observations taken at La Palma on the 5th 
of January. These are obtained from the radial velocities re¬ 
turned by HARPS pipeline once the kinematical corrections 
described above, and shown in the bottom panel of the fig¬ 
ure, are applied. 

The values do not show clear discontinuities in connec¬ 
tion with the Earth transit and suggest a complex behaviour. 
It must be noted that there is a known offset in absolute ra¬ 
dial velocities which originates from the use of the FTS solar 
spectrum as a template. This was measured in 102.5 m s~^ 
I Molaro et al.|2013 1 in coincidence of the Venus transit with 
an uncertainty of the order of few m which depends 
from the solar activity of that day. 

Both spectrographs benefit of a second fiber which sup¬ 
plies ThAr spectra simultaneous with observations and that 
can be used to correct for instrumental radial velocity drifts 
occurring over the night. The radial velocity differences with 
respect to the previous calibration provide the instrumental 
drifts for both spectrographs. 


^ Solar System Dynamics Group, Horizons Web Ephemerides 
Systems, JPL, Pasadena, GA 91109 USA http://ssd.jpl.nasa.gov 


3.1 The Radial Velocity Anomaly 

The whole set of corrected solar radial velocities obtained 
from the Jupiter moon’s spectra taken in the course of the 
3 nights from both sites is shown in Fig. [Rafter subtraction 
of the radial velocity baseline. At the beginning of the ob¬ 
servations the RV is of 107 m s~^ while at the end it is at 
108 m s~^, and we adopt here a baseline of 107.5 m for 
simplicity. The observations taken at La Palma show a sud¬ 
den drop by about 7 m s~^ after about one hour. Moreover, 
at the start of La Silla sequence the radial velocities were 
slightly lower with a difference of about 4 m between 
the two spectrographs. 

In the following day La Palma observations started at 
about mid-transit with radial velocities rising very quickly 
till they reached a peak of 37 m s“^. After the peak the 
radial velocities declined monotonically showing a break in 
the slope in correspondence to the end of the transit. The 
vertical lines in the figure mark the start, middle and end of 
the transit for Europa. To note that the peak of the radial 
velocity is reached at MJD 56662.5 in correspondence of 3/4 
of the Earth passage in the receding solar hemisphere and 
the change in the slope in declining which corresponds to 
the end of the transit. Both of them will be discussed in 
the next section where we provide an interpretation of the 
phenomenon. 

In the night following the transit we made observations 
only from La Silla. The radial velocities are back to the val¬ 
ues of the night preceding the transit. The observed pattern 
is completely at odd with our expectations. In the fraction of 
transit covered by observations the Earth was eclipsing the 
receding solar hemisphere and the RM effect should have 
produced a small blue shift of the lines as a result of the 
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Figure 5. Residuals of the difference in units of cm between 
the observed radial velocities and the RM model for the 2012 
Venus transit described in |Molaro et al.| ( |2013[ l but with the 5min 
oscillation filtered out. This illustrate the level of precision which 
can be achieved with HARPS at the sub m level in contrast 
with the radial velocity anomaly observed for the Earth’s Transit 
which is more than two orders of magnitude larger. 


prevalence of light coming from the approaching solar hemi¬ 
sphere. On the contrary we observed a change in the radial 
velocity of 37 m of opposite sign, i.e. more than two 
orders of magnitudes greater than expected. Moreover, the 
radial velocities did not show any sharp change in correspon¬ 
dence of the end of the transit. When the observations from 
the two spectrographs overlap in time the radial velocity be¬ 
haviour is similar in HARPS and HARPS-N, although there 
is a non negligible offset between the two measurements. 

The anomaly in radial velocity cannot have an instru¬ 
mental origin. This is demonstrated by the fact that the two 
observatories are giving consistent results and similar RV 
anomalies have never been observed with HARPS . An ex¬ 
ample of the precision which can be achieved in radial veloc¬ 
ities with HARPS are the observations of the Venus transit 
of 2012, which were taken with the same technique adopted 
here. For the Venus transit we obtained a remarkable agree¬ 
ment between the predicted RM effect and observations. In 
Fig. [5] the difference between the RM model computed for 
the Venus passage described in Molaro et al. (20131 and the 
observations are plotted after the observations were filtered 
for the 5m solar oscillations. The residuals of the observa¬ 
tions versus the model are of 0.55 m while the absolute 
difference is of -2 cm s~^, a difference which is within one 
sigma of the error in the normalization of the observations 
with the radial velocities observed after the transit. These 
observations were treated in the same way as those we are 
dealing here showing that large anomalies in radial veloci¬ 
ties from HARPS are not plausible. Moreover, inspection of 
asteroid observations taken with HARPS in its life span of 
12 years show that RV deviates from the mean by no more 
than ~ 5 m s~^ . Such deviations are likely correlated with 
the solar magnetic activity as can be inferred from the pres¬ 


ence of solar spots and plages on the solar surface (Lanza 
et al.|2015| ). 


Solar spots could also affect the radial velocity of the 


solar lines and indeed in Fig the solar image of 5 January 
shows the presence of several solar spots which could con¬ 
tribute at the level of few m The characteristic change 
is on a time scale of solar rotation and no effect is expected 
during the relatively short duration of the Earth transit. The 
radial velocity baseline before and after the transit also in¬ 
cludes any contributions originated by the presence of these 
solar spots. To check if short-time strong solar activity oc¬ 
curred in coincidence of the transit we inspected the Birm¬ 
ingham Solar Oscillations Network (BiSON) archives con¬ 
taining solar velocity residuals in the first days of January 
2014. The data were captured from the sites of Narrabri, 
New South Wales, Australia, Carnarvon, Western Australia, 
Izana, Tenerife and Las Campanas, Chile and provide a con¬ 
tinues monitoring of the solar activity in proximity of the 
event. The other two sites of Los Angeles and South Africa, 
were offline in those days due to bad weather. The BiSON 
velocity residuals in Fig. do not show any anomaly at the 
level observed, and suggest that the anomaly in RV we de¬ 
tected does not depend from an anomalous activity of the 
Sun. In the next sections we will see that according to our 
proposed explanation it is not a surprise that BiSON does 
not see the radial velocity anomaly. 

The effects of a microlensing onto the RM effect in the 
case of transiting planets has been studied in detail by |Os-| 


hagh et al. (20131. The RM can vanish in the extreme cases 


of particularly massive planets, but it has never been found 
to be inverted as we observed. Moreover, the size of the Ein¬ 
stein ring due to Earth observed from Jupiter is of only 47 
km which is not expected to produce any significant atten¬ 
uation of the RM effect. 

Therefore, we think that this effect is real, and we sug¬ 
gest it is due to the Opposition Surge onto the icy Europa 
as we argue in detail in the next sections after a brief intro¬ 
duction to the nature of this effect. 


3.2 The Opposition Surge effect 


The Opposition Surge is a brightening of a rocky celestial 
surface when it is observed at opposition. The increase in 
brightness is a function of phase angle and gets greater and 
greater as its phase angle of observation (f) approaches zero. 
The existence of the opposition surge was first recorded by 


Gehrels (19561 but the precise physical origin is not yet com¬ 


pletely understood and shadow hiding or coherent backscat- 
ter have been proposed. 

The former stems from the fact that when the light hits 
a rough surface at a small phase angles all shadows decrease 
and the object is illuminated by its largest extent. It was 
Hugo von Seeliger who back in 1887 explained the increase 
in albedo of Saturn’s rings to the corresponding reduction of 
the shadows on the dust particles of the rings at opposition. 

In the coherent backscatter theory the increase in 
brightness is due to a constructive combination of the light 
reflected from the surface and by dust particles. The con¬ 
structive combination is achieved when the size of the scat- 
terers in the surface of the body is comparable to the wave¬ 
length of light. At zero phase the light paths will construc¬ 
tively interfere resulting in an increase of the intensity while 
as the phase angle increases the constructive interference 
decreases. Coherent backscatter has been observed in radio 
wavelengths and detailed physical models are presented in 
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Figure 6. BiSON solar observations archival data containing ve¬ 
locity residuals in the same days of our observations. The data are 
from Narrabri, New South Wales, Australia (red points), Carnar¬ 
von, Western Australia, (black), Izana, Tenerife (red) and Las 
Campanas, Chile (green). Courtesy of Steven Hale. 


Hapke, Nelson & Smythe (19931; Hapke (2002 I and Shkura- 
tov k. Helfenstein| ( |2001^ 


It is also possible that both coherent backscatter and 
shadow hiding are operating. Which mechanism is dominant 
depends on the physical properties of the surface such as 
porosity, the mean free path, and the single particle albedo. 
Currently, theory is unable to predict the amplitudes for ei¬ 
ther mechanism (Schaefer, Rabinowitz & Tourtellotte 20091. 
Considering both explanations the Opposition Surge is also 
known as the Seeliger-Hapke effect. 


3.3 An Inverse Rossitter-McLaughlin effect 

In the following we argue that the Opposition Surge can ex¬ 
plain the radial velocity anomaly observed in proximity of 
the Earth Transit. A characteristic feature of the Opposition 
Surge is the brightening of the planet as the phase angle (j> 
decreases. Solar photons which graze the Earth have smaller 
angles than photons coming from regions of the solar disk far 
away from the Earth edge. Thus they produce an effective 
increase in the radiation coming from the region of the Sun 
just behind the Earth as it moves across the face of the Sun. 
Along its passage the Earth acts as a lens and the light mag¬ 
nification produces a radial velocity drift which is opposite 
in sign to that expected from a Rossiter-McLaughling effect, 
but of identical physical origin. The enhancement of a por¬ 
tion of the solar disk produces a distortion in the solar line 
profiles with an asymmetric contribution from the two solar 
hemispheres of the same kind of the Rossiter-McLaughlin. 
The opposite sign is because instead of an occultation there 
is an enhancement of the emission in a restricted area of 
the solar surface. Instead of receiving less radiation from 
the hemisphere the Earth is crossing, due to its occulta¬ 
tion of the solar disk we are receiving more radiation from 


it because of the enhancement produced by the Opposition 
Surge effect of the reflecting body. This effect not only com¬ 
pensates the effect of the partial solar eclipse by Earth but 
is able to produce an opposite radial velocity drift by orders 
of magnitude stronger. 

Opposition Surge has been observed in Jupiter’s moon 
Europa and has become prominent for phase angles less 
than 0 < 1 degree ( Simonelli fc Buratti||2004 1. The Jovian 
moon has a comparatively young surface rich in water ice 
which produces a high albedo. In these conditions coher¬ 
ent backscatter is expected to dominate over shadow hid¬ 
ing. However, near infrared Cassini observations have been 
interpreted as the Opposition Surge cannot be produced by 
coherent backscatter alone, but that it must have a signifi¬ 
cant shadow hiding component even in the presence of high 
albedo ( Simonelli fc Buratti|20d4 l. 

The Opposition Surge is not fully understood and we 
cannot make a quantitative prediction of the distribution of 
the light enhancement as a function of the angular distance 
from the Earth position. 

However, a simplified model which accounts for the 
asymmetrical emission from the two rotating solar hemi¬ 
spheres can explain most of the features of the RV curve we 
observed. 

We considered an area around the Earth with uniform 
enhanced emission and we computed the effect in RV as if 
it were due to the RM effect. The sign of the radial ve¬ 
locity drift is reversed to simulate the emission instead of 
the eclipse. The theoretical radial velocity anomaly of the 
Sun during the transit is computed using the formalism of 


Gimenez (20061. Since there is a degeneracy between the 


area and the intensity of emission, we just scaled the ra¬ 
dial velocity to the observed one but preserved the shape. 
The scaling factor provides the amount of light enhancement 
which is necessary for a given area assuming a uniform emis¬ 
sion, while it is very likely that it changes within the area 
as a function of the phase angle. 

We assumed the rotational velocity of the Sun Vrot is: 


uj = a+ h sin^ <;!> -I- c sin'* 0 


(3) 


where ui is the solar angular velocity measured in °/day, (j) is 
the solar latitude, and a, b, c are the coefficients derived from 
the magnetic field pattern (a = 14.37, b = —2.3, c = —1.62). 
The corresponding rotational velocity at latitude (j) defined 
by the Earth trajectory is: 


Vrot = 2-kRq ■ {a + fesin^ 4> + csin* <j>) 


(4) 


The limb darkening coefficients of the Sun are Ua = 
0.5524 and wj = 0.3637, taken from the tables of |Claret| 
(20041, for the g filter and an Atlas model for the Sun with 


solar metallicity, Tefif = 5750 K, logg = 4.5, ^ = 1 km s~ . 

The theoretical variation of the solar radial velocity dur¬ 
ing the transit computed with the above derived parameters 
is plotted as a thin line in Fig.[^ In the figure it is possible to 
see that the radial velocity anomaly does not end abruptly 
with the end of the Earth transit, but it extends further 
after it. This is not surprising since the Opposition Surge 
does not last for the time of eclipsing transit but it is also 
present when the Earth has just left or is approaching the 
solar disk, provided the solar rays are coming from portions 
of the solar disk which are at angles small enough to produce 
the Opposition Surge. For many hours after the end of the 
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Figure 7. Model of the inverse Rossiter-McLaughlin RV drift induced by an increase in the solar emissivity in the region behind the 
Earth trajectory due to the Opposition Surge. The thin black line shows the the drift expected from a small area moving with the Earth, 
while the thick line shows a drift coming from a larger area with radius of Ri 10 arcminutes required to match the start of the radial 
velocity drift. The gap between the two models is due to the numerical impossibility to compute a RM effect for total eclipse. The data 
points from la Silla (blue points) are delayed by 0.1468 MJD to compensate for the longitude difference between the observatories of La 
Palma and La Silla. Cfr text for the details. The blue dotted line shows how the expected Rossiter-McLaughlin effect but amplified by 
a factor of 50 since the real one would have the thickness of the line. The inverse RM effect detected is about 400 times larger than the 
expected RM due to the Earth transit. 


Transit, the Opposition Surge makes the solar hemisphere 
just left by the Earth brighter than the most distant one. 
Thus, the radial velocity is decreasing smoothly while the 
Earth is moving away and the phase angle is increasing. 

We observe the phenomenon extending after the transit 
but not its end since RV is still high about six hours after 
the end of the transit. It is only on the following night that 
we measured again a constant radial velocity. For symmetry 
we can assume that the Opposition Surge should also have 
started many hours before the formal start of the Earth tran¬ 
sit in coincidence with the sudden drop in radial velocities 
by about 7 m s~^ observed on 4 January at MJD 5666.204- 
5666.206 from La Palma observations of both Europa and 
Ganymede. Thus, the Opposition Surge effect started to be 
effective and produced a radial velocity change at something 
about 15 hours before the start of the Earth transit when 
the Earth was at a projected distance of about 7 arcminutes 
from the solar edge. We emphasize that it is only the differ¬ 
ence between the contributions of light coming from the two 
solar hemispheres that matters. An Opposition Surge which 
provides equal enhancement of the two solar hemispheres 
would produce a brightening but not any detectable radial 
velocity change. For symmetry the radial velocity anomaly 
should have ended also 15 hours after the end of the transit, 
i.e. in a period which is not covered by our observations. 

In our simplified model we have considered a circular 
region centered on the Earth and a radius of 6 arcminutes, 
namely 165 times the projected radius of the Earth as seen 


from Europa. In Fig. the computed radial velocities are 
overplotted to the observations approximately covering the 
transit after scaling down the RM intensity by a factor of 30. 
The predicted radial velocity rise follows the observations 
quite well, though it is somewhat less steep. The peak is 
reached when the Earth is approximately at about 3/4 of 
the solar receding hemisphere. This is the position where 
we expect the stronger effect on radial velocity due to the 
combined effect of the almost tangential rotational velocity 
and of the limb darkening of the Sun. During the decline 
a break in the slope with a more gentle decline is observed 
in proximity of the Earth egress. The region with enhanced 
emission has been enlarged to 10 arcminutes to allow the 
RV anomaly to extend well outside the transit. 

The first half of the transit could not be observed either 
from la Palma or La Silla and the observations cannot track 
the passage of the Earth in front of the approaching solar 
hemisphere where the Opposition Surge should have pro¬ 
duced a symmetrically negative radial velocity behaviour. 
Simultaneous observations from the two observatories give 
slightly different radial velocities. Those from la Silla are 
always lower than those from la Palma (see Fig.|^. The dif¬ 
ference is of about 4 m in the first night, and of about 10 
m s~^ at the beginning of the second night, but they slightly 
decrease to few m s~^ as the event faded away. While we 
cannot exclude some systematic offset between the two tele¬ 
scopes at the level of few m s~^ the difference observed 
during the opposition seems a bit too large to be explained 
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only with this systematic. Thus, it is quite possible that the 
different locations on Earth of the two observatories do not 
see exactly the same Opposition Surge. In particular, the 
distance from the Earth edge could have been relevant in 
determining the Opposition Surge intensity and therefore 
the radial velocity value. The difference between the longi¬ 
tudes of La Palma ( 28° 42.89’ N, 17° 54.29’ W ) and La 
Silla (-29 15.67’ S, 70° 43.88’ W ) is of 0.1468 MJD, while 
the distances from the equator is very similar. This means 
that after a time of 0.1468 MJD La Silla will be at approxi¬ 
mately the same distance from the Earth edge as La Palma. 
In Fig. we have shifted the data points from La Silla by 
this time difference and they provide a much better continu¬ 
ity and overlap with the values measured at La Palma. To 
note that this is achieved regardless of the fact that align¬ 
ment of the Earth and of the Jovian systems has slightly 
changed in the meantime. This would imply that the inten¬ 
sity of the Opposition Surge is very sensitive to the phase 
angle and therefore to the location of the observer on Earth, 
in particular to its distance to the Earth’s projected edges. 

It is interesting to note the possible presence of a double 
peak in proximity of the maximum of the radial velocity, 
which is suggestive of the presence of two components. While 
the broad one could be associated with a diffuse area of 
enhanced emission as we have discussed above, the latter 
narrower one could be due to a stronger emission located 
in proximity of the Earth. The result of an emission from a 
relatively small area in proximity of the Earth is plotted in 
Fig. [7] as a thin line which reproduces quite well the peak 
with a small delay of -I- 0.01 MJD. 

As we noted above no radial velocity anomalies were 
observed during the Venus Transit of 6 June 2012. The Moon 
was in opposition at about 8 degrees ahead of the Earth at 
a phase angle large enough to avoid the Opposition Surge. 


Yokota et al. (19991; Buratti, Hillier & Wang (19961 with 


their study of Clementine data estimated a 30-40% increase 
in brightness of the Moon when going from 4 to 0 degrees 
of phase angle. However, even if present, this should not 
have produced a radial velocity anomaly since it would not 
have been connected to the Venus passage in front of the 
Sun. For similar reasons the BiSON measurements obtained 
from a direct watch of the Sun do not see the radial velocity 
anomaly which is produced by the magnification of a portion 
of the Sun induced by the Opposition Surge on the Earth 
passage. 

We note also that the presence of a strong Opposition 
Surge during the Earth transit is probably the explanation 
of the lack of detection of the luminosity drop in the flux 
due to the Earth occultation which has been searched for 
unsuccessfully by many teams. 


4 CONCLUSIONS 


We followed the Earth transit of 5 Jan 2014 as seen from 
Jupiter by means of observations of Jupiter’s moons Europa 
and Ganymede. The observations were made with HARPS 
spectrograph at La Silla, Chile, and with HARPS-N spec¬ 
trographs from La Palma, Canary Islands, originally aimed 
to detect the RM effect due to the Earth passage on the 
face of the Sun. We followed the same technique success¬ 
fully adopted for the 2012 Venus Transit (Molaro et al. 20131 


where the RM effect was measured and found in agreement 
with the theoretical model within few cm s“^. In the case of 
the Earth transit the expected modulation in radial veloc¬ 
ities was of K, 20 cm s“^. Instead, an anomalous and very 
large radial velocity drift was observed. The half amplitude 
of the radial velocity drift observed was as large as 35 m s“^, 
i.e. about four hundred times higher and opposite in sign. 

The similar behaviour in the observations taken from 
both telescopes rules out an instrumental origin and sug¬ 
gests a physical origin which we identified as the product of 
the Opposition Surge effect onto Europa’s icy surface. The 
Opposition Surge effect amplifies the intensity of the solar 
radiation from the portion of the Sun crossed by the Earth 
and produces a sort of inverse Rossiter-McLaughling. This 
phenomenon has never been observed before and is associ¬ 
ated to the rather unique geometry in which we observed 
the Earth transit. In fact, simultaneous radial velocities ob¬ 
tained by BiSON through a direct solar watch do not show 
the radial velocity anomaly and rule out that they originate 
in the Sun. 

A toy model which assumes an enhancement of the so¬ 
lar radiation from a projected solar region centered on the 
Earth’s position produced by the Opposition Surge explains 
the general behaviour shown by the radial velocity measure¬ 
ments. In particular, we are able to explain why the anomaly 
is observed also before and after the Earth Transit, and the 
differences in radial velocities measured by the two observa¬ 
tories as due to the different distances from the Earth edge, 
as well as the presence of a second peak associated with the 
smaller projected solar region around the Earth but with 
greater intensity, and why we did not see a similar anomaly 
in the 2012 Venus Transit. 

The Opposition Surge effect provides a coherent and 
plausible description of the anomaly in radial velocity as an 
inverse Rossitter-McLaughlin that we observed for the first 
time during the Earth Transit. The effect could be observed 
again every time the Earth is seen in transit against the Sun 
from other planets or smaller bodies in the solar system. 

The next Earth transit will occur from jupiter in 2026, 
but it will be a grazing transit quite unfavourable to any 
kind of observations ( Meeus|[l989 |. However, since we have 
observed the effect of the Opposition Surge when the Earth 
was at an angle as high as about 10 arcminutes, we can 
predict that this same phenomenon can be observed again 
although with a minor amplitude in radial velocities. 
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